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Abstract 

Linear precoding is a relatively simple method of MIMO signaling that can also be optimal in 
certain special cases. This paper is dedicated to high-SNR analysis of MIMO linear precoding. The 
Diversity-Multiplexing Tradeoff (DMT) of a number of linear precoders is analyzed. Furthermore, since 
the diversity at finite rate (also known as the fixed-rate regime, corresponding to multiplexing gain of 
zero) does not always follow from the DMT, linear precoders are also analyzed for their diversity at fixed 
rates. In several cases, the diversity at multiplexing gain of zero is found not to be unique, but rather 
to depend on spectral efficiency. The analysis includes the zero-forcing (ZF), regularized ZF, matched 
filtering and Wiener filtering precoders. We calculate the DMT of ZF precoding under two common 
design approaches, namely maximizing the throughput and minimizing the transmit power. It is shown 
that regularized ZF (RZF) or Matched filter (MF) suffer from error floors for all positive multiplexing 
gains. However, in the fixed rate regime, RZF and MF precoding achieve full diversity up to a certain 
spectral efficiency and zero diversity at rates above it. When the regularization parameter in the RZF is 
optimized in the MMSE sense, the structure is known as the Wiener precoder which in the fixed-rate 
regime is shown to have diversity that depends not only on the number of antennas, but also on the 
spectral efficiency. The diversity in the presence of both precoding and equalization is also analyzed. 

I. Introduction 

Precoding is a preprocessing technique that exploits channel-state information at the transmitter (CSIT) 
to match the transmission to the instantaneous channel conditions [[T], 0, |[3l- Linear and non-linear 
precoding designs are available in the literature |4|. Linear precoding in particular provides a simple and 
efficient method to utilize CSIT. Linear precoding has been shown to be optimal in certain situations 
involving partial CSIT Q, Q, however, in many instances the main motivation of linear precoders is to 
simplify the MIMO receiver. 
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USA, E-mail: ahmed.mehana@student.utdallas.edu; aria@utdallas.edu. 
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Lineal- precoders include zero-forcing (ZF), matched filtering (MF), Wiener filtering, and regularized 
zero-forcing (RZF). The ZF precoding schemes were extensively studied in multiuser systems as the ZF 
decouples the multiuser channel into independent single-user channels and has been shown to achieve a 
large portion of dirty paper coding capacity [V]. ZF precoding often involves channel inversion, using 
the pseudo-inverse of the channel or other generalized inverses f?]. Matched filter (MF) precoding [8J, 
similarly to the MF receiver, is interference limited at high SNR but it outperforms the ZF precoder at 
low SNR 14). The regularized ZF precoder, as the name implies, introduces a regularization parameter 
in channel inversion. If the regularization parameter is inversely proportional to SNR, the RZF of |[9l is 
identical to the Wiener filter precoding lH. Peel et al. Q introduce a vector perturbation technique to 
reduce the transmit power of the RZF method, showing that in this way RZF can operate near channel 
capacity. 

This paper analyzes the diversity of MIMO linear precoding with or without linear receivers. We show 
that a MIMO ZF precoder with a maximum likelihood receiver has minimal spatial diversity, and that 
Wiener precoders produce a diversity that is a complex function of spectral efficiency and the number 
of transmit and receive antennas. At very low rates, the Wiener precoder enjoys a maximal diversity 
which is the product of the number of transmit and receive antennas, while at very high rates it achieves 
a minimal diversity which is the same as ZF diversity. These results are reminiscent of MIMO linear 
equalizers |[TOl . even though in general the behavior of equalizers (receive side) can be very different 
from precoders (transmit side) and the analysis does not carry from one to the other. We also show that 
MIMO systems with RZF and MF precoders (together with optimal receivers) exhibit a new kind of 
rate-dependent diversity that has not to date been observed or reported, i.e., they either have full diversity 
or zero diversity (error floor) depending on the operating spectral efficiency R. 

We also provide DMT analysis for all precoders mentioned above. The fact that DMT and the diversity 
under fixed-rate regime require separate analyses has been established for MIMO linear equalizers [11], 
lITOl and is by now a well-understood phenomenon. Essentially, the reason is that various fixed rates 
(spectral efficiencies) for MIMO precoding result in distinctly different diversities, whereas DMT analysis 
assigns only a single value of diversity to all fixed rates (all fixed rates correspond to multiplexing gain 
zero). 

Remark 1: Due to symbolic similarities, it may be tempting to draw the conclusion that if d{r) 
is the diversity at multiplexing gain r, then substituting r = in the same mathematical expression 
will give the diversity at multiplexing gain zero d{Q). However, despite appearances, there is no solid 
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relationship between d{r) and d{0). The standard DMT arguments are based on the seminal work of 
Zheng and Tse |[T2l whose developments depend critically on the positivity of r. For example, the proof 
of lllli Lemma 5] depends critically on r being strictly positive.More importantly, the asymptotic outage 
calculations in |[T2l p. 1079] implicitly use r > and result in the outage region: 

A = {a: ^(1 - Qi)+ < r} 

i 

where are the exponential order of the channel eigenvalues, i.e., Aj = p~°'\ If we set r = this 
expression implies that the outage region is always empty, which is clearly not true. 

Thus, the DMT as calculated by the standard methods of |12| does not extend to r = 0. The DMT 
d{r) is sometimes continuous at zero, including e.g. the examples in |12|, but continuity at r = does 
not always hold. In fact, there are systems where d(0), the diversity at multiplexing gain zero, is not even 
uniquely defined. It is possible for diversity to take multiple values as a function of rate R. This fact 
has been observed and analyzed, e.g., in lIlTI . fTOl . |fT31 . The work in the present paper also produces 
several examples of this phenomenon. 

This paper is organized as follows. Section |ll] describes the system model. Section HID provides outage 
analysis of many precoded MIMO systems. Section |IV] provides the DMT analysis. The case of joint 
linear transmit and receive filters is discussed in Section IVl Section IVllprovides simulations that illuminate 
our results. 

II. System Model 

A MIMO system with linear precoding is depicted in Fig. [T] This system uses the linear precoder 
to manage the interference between the streams in a MIMO system to avoid a lattice decoder in the 
receiver. We consider a fiat fading channel H G C^^^^, where M and N are the number of transmit 
and receive antennas, respectively. While M ^ N when using linear precoding alone, we have N ^ M 
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or M ^ N when using preceding together with receive-side linear equalization depending on whether 
the precoder is designed for the equalized channel or the equalizer is designed for the precoded channel 
(see Figure O. The input-output system model for flat fading MIMO channel with M transmit and 
receive antennas is given by 

y = WHTx + Wn (1) 

where T G c^^^^ is the precoder matrix, W G Q^xN jj^g receiver side equalizer. The latter may be 
set to identity in cases where the receiver does not use linear equalization. The number of information 
symbols is B ^ min(M, N), the transmitted vector is x G C and n G C is the Gaussian noise 
vector. The vectors x and n are assumed independent. 

We aim to characterize the diversity gain, d{R, M, N), as a function of the spectral efficiency R 
(bits/sec/Hz) and the number of transmit and receive antennas. This requires a Pairwise Error Probability 
(PEP) analysis which is not directly tractable. Instead, we find the exponential order of outage probability 
and then demonstrate that outage and PEP exhibit identical exponential orders. 

The objective of linear preceding (possibly together with linear equalization at the receiver) is to 
transform the MIMO channel into min(M, N) parallel channels that can be described by 

yk = V^Xk + rik, k = l,...,B (2) 

where 7^ is the SINR at the k-th receiver output and B = min{M, N). Following the notation of fT4ll . 
we define the outage-type quantities 

Po«t (i?, iV, M) ^ P(/(x; y ) < R) (3) 

dout{R,N,M) = - hm (4) 

p^oo log p 

where p is the transmitted equivalent SNR. 

The outage probabilities of MIMO systems under joint spatial encoding is respectively given by ifTTI . 
HI 



P ^ 



J]log(l + 7fc) ^ii] (5) 
k=i ^ 



We shall perform outage analysis for different precoders/equalizers as the first step towards deriving the 
diversity function. We then provide lower and upper bounds on error probability via outage probabilities. 
This two-step approach was first proposed in |[T2l due to the intractability of the direct PEP analysis for 
many MIMO architectures. 
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We denote the exponential equality of two functions f{p) and g{p) as f{p) = g{p) when 

lim = lim 

p-s>oo log(p) p^oo log(p) 

In the following, we shall need to specify various upper and lower bounds or approximations of the 
SINR 7, which will give rise to a number of pseudo-SINR variables 7, 7, and 7. 

III. Precoding Diversity 

In this section we analyze a linearly precoded MIMO system where M > N and the number of data 
streams B is equal to N. 

A. Zero-Forcing Precoding 

The ZF precoder completely eliminates the interference at the receiver. ZF precoding is well studied 
in the literature via performance measures such as throughput and fairness under a total (or per antenna) 
power constraint HH and references therein]. 

1 ) Design Method I: One approach to design the ZF precoder is to solve the following problem H 

T =argrmnE[||Tx||^] (6) 
subject to HT = I 

The resulting ZF transmit filter is given by 

T = /? H^(HH^)-i G C^'^^^ (7) 

where /3 is a scaling factor to satisfy the transmit power constraint, that is (4] 

/3hv{TT^) ^ p (8) 

where we assume that the noise power is one and that the information streams are independent. From ([D, 
the received SINR per stream is thus given by 



^^FP = P (9) 

tr(HH^^^-i ■ ^ ' 



Using the outage probabiUty is given by 

tr((HH 



Pout = P(iVlog(l + rT77T^Tf7wTT) ^^1 (10) 
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A direct evaluation of (ITOl) is intractable since the diagonal elements of (HH^)^^ are distributed 
according to the inverse-chi-square distribution |[T6l . ifTTIl . We instead bound (fTOl) from below and above 
and show that the two bounds match asymptotically. 

Let Afc be the k-th eigenvalue of HH^. Equation ( fTOl ) can be written as 



Z-^k — 1 A/s 



which can be bounded as 



nut^P(^A^log(l + ^A„,in)^i?j (11) 

= p(^A^in^iV(2#-l)i2p-ij 

= P(A„iin ^P"'). (12) 

The marginal distribution /i(A) of Amin is /i(A) = cA^*'^"^) ifTTl where c is a constant, therefore the 
bound in (fT2l) can be evaluated yielding: 

Pout ^ (13) 

We now proceed with a lower bound on outage. The outage probabiUty in (ITOl ) can be bounded: 



^ Pf iVlog(l + '^7^) ^ R\ 



fcfc 

z ^ p-') (14) 

where we have made a change of variable z = . 

The random variable z in ([141) is distributed according to the chi-square distribution with 2(M — A^+1) 
degree of freedom, i.e. z ~ X'^2{M — N + 1) ||T6ll . Thus the bound in ([T4l) can be evaluated [11] yielding: 

Pout^ (15) 

From ([T3] ) and (fTSl ). we conclude that the diversity of MIMO system using the ZF precoder given 
by ^ and joint spatial encoding is 

d^PP = M -N + 1. (16) 
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B. Zero-Forcing Precoding: Design Method II 

Notice that the ZF precoder design in ^ minimizes the transmitted power. Another approach for 
ZF precoding design allocates unequal power levels across the transmit antennas to optimize some 
performance measure. For instance, consider the optimization problem ifTSl 

max f{pk) 

subject to HT = diag{^, ...,^^} 

E||Tx|p^p (17) 

where f{pk) is an arbitrary function of the transmitted power pk on the k-th antenna. 

The optimal solution for (fTTl ) (assuming independent transmit signaling) is given by |T5^, Theorem 1]: 

T = H^(HH^)^idiag{v^,...,Vm} (18) 
where pk is obtained by solving 

max f{pk) 
Pk 

subject to Pk [(HH^)-i] ^ p (19) 

k 

In our case, we maximize the throughput, therefore f{pk) = X]fclos(l + ^k^^)- After setting the 
derivatives of the appropriate Lagrangian function to zero, the solution of the power allocation problem 
in ( fT9l ) is given by 

„ _ P + Efc=i(HH^)fcfc 

" — jim%i — ' ™ 

Substituting (l20l ) in ([S]), the outage probability is given by 

N . 



N\ogY,{^^^-^)^-R] (23) 



k=l p 



N 

k=l 



" 1 

Y — 
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p-iM-N+i)_ (25) 



where the exponential equaUty (|2Tl) holds at high SNR, (1231) follows from Jensen's inequality, and the 
transition from (l24l) to (l25l ) again due to the marginal distribution of Amin via the method of llT3l . 

A lower bound on the outage probability can be given as follows. Starting with ((2T]) and using Jensen's 
inequality we have 

N 

P 



Pont= P Vlogf ^^^) ^ R 



k=l -'V-"- Jkk 

N 



>P(iVlog-^V ^^^^^V (26) 



The singular value decomposition of H and the corresponding eigen decomposition of HH^ are given 

by 

H =urv^ 

HH^ = UAU^ 

where U G £_NxN y g ([^MxM unitary matrices, T G M^^*^ is a rectangular matrix with non- 
negative real diagonal elements and zero off-diagonal elements, and A = FT^ G M^^^ is a diagonal 
matrix whose diagonal elements are the eigenvalues of HH^. Let be the A;-th column of U^. We 
have 



N 



|2 



1=1 

where Uki is the {k, I) entry of the matrix U. 
The bound in (l26l ) can be rewritten 

N 



{UU^rJ = u^A-W = ^ (27) 



^ I— 1 2^1=1 ^ 



k=\ L^l=\ p\i 

N 



k=l 1^1=1 1+pAi 



j\r2 \uui\^- 

k=i 2^1=1 i+px 

We can lower bound the probability in (1281) in by observing that the term 

N 



ly L 



- - 1=1 l+p\i 
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is similar to |[T3l Eq.(18)], thus the analysis of |13| appUes and we obtain 

Po.t^p(An.i„^p-l) =p-iM-N+l)_ (29) 

Thus, the MIMO ZF precoding with unequal power allocation ( [T9l ) achieves diversity order M — + 



Recall that the diversity is defined based on the error probability. In Appendix|A]we provide the pairwise 
error probability (PEP) analysis for the zero-forcing and regularized zero-forcing precoded systems and 
show that the outage and error probabilities exhibit same diversity. 

C. Regularized Zero-Forcing Precoding 

In general, direct channel inversion performs poorly due to the singular value spread of the channel 
matrix |9]. One technique often used is to regularize the channel inversion: 

T = /3H^(HH^ + cI)^^ (30) 

where /3 is a normalization factor and c is a fixed constant. 
Recall 

y = HTx + n = /?UA(A + cI)"^U^x + n (31) 

allowing us to decompose the received waveform at each antenna into signal, interference, and noise 
terms: 

N 

Vk 



where the scaling factor /3 is given by /3 = ^ and 



/3 X] i^J^^^'^kiu*ii]xi + nk (32) 

i=l,i=^k 1=1 



r] = tr[(HH^ + cI)"^HH^(HH^ + cl)"^] 
= tr[(UAU^ + cI)"^UAU^(UAU^ + cl)"^] 
= tr[U(A + cI)-iA(A + cI)-^U^] 

= tr[A(A + cI)-2]=j;^^-^. (33) 



The received signal power is given by 



Pt = £;||HTx|p 
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E 



E 



fi'hxi UA(A + cI)"iU^xx^U(A + cI)"^AU 



phri A(A + cI)-iU^xx^U(A + cI)-UU^U 



1 ATTi^l 



/3^tr( A(A + cI)"^U^^(xx^)U(A + ciy'A 



^-tr[(A + cI)-2A2^ 
where we have used E(xx^) = 



N 



■g(A.+c) 



(34) 



The SINR is evaluated by computing the signal and interference powers from ([32]). For a given channel 
H, the power of desired and interference signals at the A;-th receive antenna are respectively given by 



p(fc) _ l^'^pf A; ,2 



,2. A, 



p(fc) ^ ^ V- V- 
^ iV . ^ \i + c 



UklUii 



i=i,ijkk 1=1 

Thus the SINR for the A;-th signal stream is given by 



(35) 
(36) 



(fc) 



7fc 



D 



V^A^ A, * 



(37) 



(38) 



recall -q is given by (|33]) . 

Defining the exponential order of eigenvalues A/ = p^"' in a manner similar to {\2 

2 



Ik 



E 



+ N p-^ri 



E 



(39) 



+ A^p-^E^IiP 



where the asymptotic equality follows because in all terms c dominates p a fact that also implies 

^ = E/P""'- 
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Multiplying the numerator and denominator of ( 1391 ) by p , we have 

2 



E 



E7V * 
/=1 UklUiiP 



(40) 



+ A^E^Ii/^^- 



■cti 



The sum in the numerator of (1401 ) is, in the SNR exponent, equivalent to: 



^l-Omi„ 



(41) 



where we use the fact that \uki\'^ = 1. Similarly, for the first term in the denominator of (l40l ) 



N 



1=1 



„2-2o„ 



i^fc Z=l 



Af 2 



Y^Wki 

i^k 



(42) 



where we define 



1=1 ^klU^l 



. Notice that Wki < 1- 



Using dHJ and (02)), the SINR in (00)) is given by 



Ik-- 



1-Q„ 



If all a£ > 1 then the exponents of p are negative and the denominator is dominated by its second 
term, which also dominates the numerator. If at least one of the < 1, then the maximum exponent 
which corresponds to amin dominates each summation. Thus we have: 



(43) 



lk = { 



a; > 1 
otherwise 



We now concentrate on the case where there exists at least one <1. We define 



k,i 
k^i 



therefore in this special case we have: 



(iV-l)(pi-'^-")'//^in + Arpi- 



(44) 



(45) 



(46) 
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A - 

= 7 



Thus in general 



A - 

= 7 



where v is a new random variable defined as: 

Ka if Ofc > 1 VA; 



1 otherwise 



where k„ = 



We can now bound the outage probability as follows 

N 

Pout 



fc=i ^ 

N . 

5]log(l + 7) ^i?) 

fe=i ^ 



(iV - 

V 



< 2^/^ - 1 



where G = (2^/^-1) (TV -1). 

The bound in dSOl ) can be evaluated as follows 



(47) 



[N - 1)/Xmin 



(49) 



^ e (50) 



^ G = P( ^ Q\v = Ko)^{v = Ka) + P( ^ G|i^ = = 1) 



A'min / A'min /^min 



= P(k„ ^ G ^„,i„)P(i/ = K„) + ^ Q)P(i/ = 1) . (51) 

/^min 

Notice that P(ko ^ G /imin) = 1 since Ka is vanishing at high SNR and G and fimm are positives. We 
now need to compute P(z^ = Ka) and P(i/ = l), or equivalently P({aA; > 1 VA;}) and its complement. 
We quote one of the results of [10]. 

Lemma 1: Let {A„} denotes the eigenvalues of a Wishart matrix HH^, where H is an x M matrix 
with i.i.d Gaussian entries, and let a„ = If la„ denotes the number of a„ that are greater than 
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one, then for any integer s ^ we have ifTOl Section III- A] Q 

P(l,^ =s) =p-(-^+{A'^-^W. (52) 

Thus setting s = N (i.e. all a„ > 1) in (l52l) yields 

F{iy = Ka) = P(la„ =N) = p-^^ (53) 

F{iy = 1) = 0(1) (54) 

where 0(1) is a non-zero constant with respect to p. 

Evaluating (ISTl ) depends on the values of G which is always real and positive. If < 1 then we have 



because ^{jf— ^0) = as l/pmm > 1. On the other hand if > 1 then 

^ ^0") =p-A^^ + p(^ ^0)O(1) (56) 

= Oil) (57) 

since ^ 0) is not a function of p because p is independent p. For the set of rates where > 1, 
equation (157] ) implies that the outage probability in (|82l ) is not function of p and thus the diversity is 
zero, i.e. the system will have error floor. The set of rates for which > 1 are 

R > Nlog{j^) ^Rth- (58) 

This concludes the calculation of a lower bound on the outage probability. A similar approach will 
yield a corresponding upper bound, as follows. Let 



/^max = max|ufc;/<;,| (59) 



A lower bound on the SINR is given as 



Ik ^ TTp— T (60) 

A ~ 

= 7- 



'Note that flOl analyzes linear MIMO receiver where it is assumed A'^ ^ M. It can be easily shown that the above Lemma[T] 
applies for the case considered here where M A'^. 
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The outage probability is bounded as 



k=l 




/^max 




(61) 



We can evaluate (|6T]) in a similar way as (ISTl) . establishing that the outage diversity = MN 

if the operating spectral efficiency R is less than Rth = N\og{j^), and (i^^ = if i? > Rth- 
This shows that the performance of RZF precoder can be much better than that of the conventional ZF 
precoder MIMO system whose diversity is M — + 1 independent of rate. 

Recall that diversity is the SNR exponent of the probability of codeword error. In Appendix |Al we 
show that the outage exponent tightly bounds the SNR exponent of the error probability. Thus we have 
the following theorem. 

Theorem 1: For an M x MIMO system that utilizes joint spatial encoding and regularized ZF 
precoder given by ( [30l ). the outage diversity is d^^^ = MN if the operating spectral efficiency R is less 



Remark 2: Rth is a monotonically decreasing function of N with the asymptotic value lim n^oo Rth = 
^^^2 ~ 1.44. Overall we have 1.44 < Rth < 2, leading to an easily remembered rule of thumb that applies 
to all antenna configurations. Regularized ZF precoders always exhibit an error floor at spectral efficiencies 
above 2 b/s/Hz, and enjoy full diversity at spectral efficiencies below 1.44 b/s/Hz. 

D. Matched FUter Precoding 

The transmit matched filter (TxMF) is introduced in HI, llH. The TxMF maximizes the signal-to- 
interference ratio (SIR) at the receiver and is optimum for high signal-to-noise-ratio scenarios llH. The 
TXMF is also proposed for non-cooperative cellular wireless network |]T8l. The TxMF is derived by 
maximizing the ratio between the power of the desired signal portion in the received signal and the 
signal power under the transmit power constraint, that is [4| 



than Rth = iVlog (j^), and d^^^ = if R > Rth ■ 




(62) 



subject to: £'||Tx|p ^ p 



where y is the noiseless received signal y = Tx. 



February 27, 2012 



DRAFT 



15 



The solution to (162]) is given by 



H 



(63) 



with 



/3 



1 



tr(H^^H) ■ 

We now analyze the diversity for the MIMO system under TxMF. The received signal is given by 

y = HH^x + n = /3UAU^x + n. 
The received signal at the A;-th antenna 



(64) 



Vk = /3( ^M^kl?']: 
^ 1=1 ' 

N , N V 

/3 ^ f ^ hukiu*ii jxi + rik 

i=lAj^k ^ 1=1 ' 



(65) 



The SINR at A;-th receive antenna is 



P AT 



7fc = - 

a2P_ 

P N l^i=l,i^k 

Substitute with the value of /3 and = 

Ik = 



{Ya=i Muki?^ 



T.i=i kukiu*i 



+ 1 



Z—/i=l,iy^k 



Elip-'^'^kiui 

Observe that (1661 ) is the same as the SINR of the RZF precoded system given by (|39l ). Hence the 
analysis in the present case follows closely that of the outage lower bound of the RZF precoder, with 
the following result: the system can achieve full diversity as long as the operating rate is less than Rt^ 
given in (l58l) . The pairwise error probability analysis is also similar to that of the RZF precoding system 
(given in Appendix |A]) which we omit for brevity. Thus we conclude that Theorem \T\ applies for the 
TxMF precoder. 



(66) 



+ Np-^EliP 



E. Wiener Filter Precoding 

The transmit Wiener filter TxWF minimizes the weighted MSE function. 

{T,/3} =argminT,^ii;(||x-r'y||2) 
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subject to ^(||Tx||2) ^ p. (67) 

Solving ( [67l ) yields 

T = /3F-^H^ (68) 



with 



F = j^H^H + — I 

where /3 can be interpreted as the optimum gain for the combined precoder and channel H. 

Notice that the TxWF precoding function is similar to that of the MMSE equalizer |[T9l . Indeed the 
SINR of both systems are equivalent. To see this, we first compute the SINR for the precoded H G £,MxN 
(with M ^ N) MIMO channel 

£||(T H)fcfc|2 

Fl(TH)fcfc|2 ^^^^ 



#Ejj(T H),,|2+tr(F-2H^^H) 
where we have used the independence of the transmitted signal to compute (ITOb . 

Now consider a MIMO channel H2 = € C^^'^'^ . The MMSE equalizer for this channel is given 

by 

We = (H^H2 + -I)-^H^. (72) 
P 

The received SINR for that system is given by 

MMSE ^ ^|(We H2)fcfcP 

#E.vJ(WeH2)HP+tr(WeWe)' 

Since We H2 = TwFpH and tr(WeWe) = tr(F^2H^^H), we conclude that 7^ ^^'^-^ = 7^^^. 
Hence the diversity analysis of (TOl, |[T3l for the MIMO MMSE receiver applies for the MIMO Wiener 
precoding system. It is shown in lITOl that this diversity is a function of rate R and number of transmit 
and receive antennas. We thus conclude the following. 

Lemma 2: Consider a channel H G £^MxN diversity of the MIMO system under Wiener filter 
precoding is given by 

^WFP ^ p^2-f ]2 + (M - iV)[iV2^#] (74) 
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where (•)+ = max(-,0) and [•]. 

Remark 3: It is commonly stated that MMSE and ZF operators "converge" at high SNR. The devel- 
opments in this paper as well as ifTTl serve to show that although not false, this comment is essentially 
fruitless because the performance of MMSE and ZF at high SNR are very different. This apparent 
incongruity is explained in the broadest sense as follows: Even though the MMSE coefficients converge to 
ZF coefficients as p — )■ oo, the high sensitivity of logarithm of errors (especially at low error probabilities) 
to coefficients is such that the convergence of MMSE to ZF coefficients is not fast enough for the logarithm 
of respective errors to converge. 

IV. Diversity-Multiplexing Tradeoff in Precoding 

For increasing sequence of SNRs, consider a corresponding sequence of codebooks C{p), designed at 
increasing rates R{p) and yielding average error probabilities Pe{p)- Then define 

r = iim 

p-s>oo log p 

p-^oo log p 

For each r the corresponding diversity d{r) is defined (with a slight abuse of notation) as the supremum 
of the diversities over all possible codebook sequences C{p). 

From the viewpoint of definitions, the traditional notion of diversity can be considered a special case 
of the DMT by setting r = 0. However, from the viewpoint of analysis, the approximations needed in 
DMT calculation make use of R{p) being a strictly increasing function, while for diversity analysis R 
is constant (not strictly increasing function of p). Thus, although sometimes DMT analysis may produce 
results that are luckily consistent with diversity analysi^j {r = 0), in other cases one may not be so 
lucky and the DMT analysis may produce results that are inconsistent with diversity analysis. Certain 
equalizers and precoders fall into the latter category. In the following, we calculate the DMT of the 
various precoders considered up to this point. 

1 ) ZF Precoding: Recall that two ZF precoding designs have been considered. For the ZF precoder 
minimizing power, given by ([7]), the outage upper bound in ([TT]) can be written as 

Pout ^ P(Amin ^ P^--^^) (75) 
^E.g. the point-to-point MIMO channel with ML decoding. 
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where we substitute R = r log p to obtain (1751) . and equation (1761 ) follows in a manner identical to the 
procedure that led to ([T3l) . 

Similarly the outage lower bound (fT4l) can be written as 

Pout > P{Z ^ p(^-l)) 

^ (A/-JV+1)(1-^) 



From (1761) and (177]) we conclude 



(M- + 1)(1 



(77) 



(78) 



The DMT of the ZF precoder maximizing the throughput, given by (ITSI ). is obtained in an essentially 
similar manner to the above, therefore the discussion is omitted in the interest of brevity. 

2) Regularized ZF Precoding: We begin by producing an outage lower bound. To do so, we start 
by the bound on the SINR of each stream k obtained in (l44b . and further bound it by discarding some 
positive terms in the denominator. 



Ik 



1 k = l 



k = 1 
k > 1 



1 



We can now bound the outage probability 



k>l 



■ i V \ 

J]log(l+7fe) ^i?) 

■ k=l ^ 

J]log(l + 7fc) ^i?] 
^ k=\ ' 

( ^ 1 

^P iVlogj;^(l+7fc)^i? 

fc=i 

^ 1 

NM, ^ ^ '"^ 



k=l 



(79) 
(80) 
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N 



■ k=l 

N 

I V 

\ l7/,;,;/7;*,.P 

k=2 



where we have used the Specht bound in ( 1791 ) in a manner similar to 1 10]. Equation (l80l ) and (ISTl ) follow 
similarly to (lOl Section III-B] 

For notational convenience define 



, A 1 -c^ 1 

= 1 n2 + / V 1 r72 



Then the bound in (l82l) can be evaluated as follows: 



= P(z^^ ^ p^|zy = 0)P(z^ = 0) +P(i/7/; ^ = l)P(z/ = 1) 

= P(0 ^ p^)P(i/ = O) + P(?/; ^ p^)P(z^ = l) 

= p-^^^ + P(^^/9«) 0(1). (83) 

^ p-A^^ + 0(1) (84) 

= 0(1) (85) 

where (1831 ) follows from Lemma [U and (l84l ) is true as long as P(V' ^ = the proof of which 

is relegated to Appendix |B] 

Since the outage lower bound (l84l ) is not a function of p, the system will always have an error floor. 
In other words the DMT is given by 

d^^^^ir) =0 0<r<B (86) 

We saw earlier that in the fixed-rate regime RZF precoding enjoys full diversity for spectral efficiencies 
below a certain threshold, but it now appears that DMT shows only zero diversity. DMT is not capable 
of predicting the complex behavior at r = because the DMT framework only assigns a single value 
diversity to all distinct spectral efficiencies at r = 0. A similar behavior was observed and analyzed for 
the MMSE MIMO receiver EH, ||T3l, ifTOl . 
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Fig. 2. MIMO with linear precoder with receive-side equalization 



3) Matched Filter Precoding: The DMT of the MIMO system with TxMF is the same as the DMT 
given by ( [86l ) due to the similarity in the outage analysis (see Section IIII-DI ). We omit the details for 
brevity. 

4) Wiener Filter Precoding: Since the the received SINR of the MIMO system using TxWF precoding 
is the same as that of MIMO MMSE receiver, we conclude from [TBI that the DMT for the TxWF 
precoding system is 

d^^^(r) = (M-iV + l)(l- ^) + . (87) 

Similarly to the MIMO MMSE receiver |[T3l, HOl, we observe that DMT for the MIMO system with 
TxWF does not always predict the diversity in the fixed rate regime given by (1741) . 

V. Equalization for Linearly Precoded Transmission 

The objective of a precoded transmitter is to separate the data streams at the receiver. In other words, 
linear precoding is a method of interference management at the transmitter. In general, precoded systems 
do not require interference management at the receiver, however, once a transmitter is designed and 
standardized (as precoders have been), some standards-compliant receivers may opt to further equalize 
the precoded channel (see Figure |2]). This section analyzes the equalization of precoded transmissions. 

When the transmit and receive filters can be designed jointly and from scratch, singular value decom- 
position becomes an attractive option whose diversity has been analyzed in 1201 . The distinction of the 
systems analyzed in this section is that the precoders can be used with or without the receive filters, 
while with the SVD solution neither the transmit nor the receive filters can operate without each other. 

A snapshot of some of the results of this section is as follows. It is shown that equalization at the 
receiver can alleviate the error floor that was observed in matched filter precoding as well as regularized 
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ZF preceding. It is shown that MMSE equaUzation does not affect the diversity of Wiener filter precoding, 
but ZF equaUzation does indeed affect the diversity of Wiener filter precoding in a negative way. 

Recall that in the system model given in Section|lI]we have defined the precoder and equalizer matrices 
T € C'^'^^^ and W € C^^^, respectively, where B is the number of information symbols not to exceed 
min(M, N). In most wireless systems, the equalizer at the receiver is designed to equalize the compound 
channel (HT) composed of the precoder and the channel (rather than designing the precoder for the 
equalized channel (WH) although it is possible). In such case we have M ^ N and we set B = N. 

A. ZF Equalizer 

The ZF equalizer is analyzed when operating together with various precoders, as follows. 
1 ) Wiener Filter Precoding: The TxWF precoder is given by 



= pn^fuii^ + —In] (88) 

P V P J 

where (l88l) follows from 111] Fact 2.16.16] q The scalar coefficient /3 is given in (l69l) and, similar to (1331) . 
it can be written as /3 = 1/ 

The ZF equalizer for the precoder and the channel is given by 

WzF = (M^Uy^m^ (89) 

The composite channel H is given by 

H = HT. 

The received signal is given by 

y = WzFiiTx + WzFn. (90) 

The filtered noise n = Wzft^ is is a complex Gaussian vector with zero-mean and covariance matrix 
Rn given by 

Rfi = [M^H]-i 

^Let A G C"""" and B G C"""" then (I„ + AB)-^A = A(I„ + BA)"\ This fact can be proved via Matrix Inversion 
Lemma. 
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= [(HH^ + Np-^ I)-^(HH^)2(HH^ + N p-^ l)-^] 

= [UA(A + Np-^ I)-iU^UA(A + Np-^ ly^U^] 

= [UA2(A + Afp-il)-2u^]"' 

where we have used the eigen decomposition HH^ = UAU^. The noise variance of the output stream 
k is therefore 

i?.(fc,fc) = ^(^^^^)V.d^ (91) 

where (|9l1 ) follows in a similar manner as ( [27l ). We can compute the signal-to-noise ratio of the ZF filter 
output: 



NRn{k,k) 

(92) 



Due to the complexity of (l92l ) we proceed to bound the outage from above and below. The upper 
bound on outage is calculated as follows. Since \uki\ ^ 1, 

> PJ^ (93) 



l/N 

EN v^Af / pi-°i+jV \ S 



(94) 



A - 



7- (95) 
where we have substituted A; = p^"' in ( |94l ). Thus the outage probability is bounded as 

N 

Pom. 

N 



^prj]iog(i+7) ^i?] 

^ k=l ' 



7 ^ 2- - Ij (96) 

Similarly to the previous analysis, we examine the SINR bound 7 for different values of a/. Define the 
set ;B = {/ I > 1} and the event 

(i = {\B\ = N} (97) 

we have 

Pout^P(7^2#-l 
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P(£) +P 7 ^ 2- - l|£ P(£) (98) 



^ P 7 ^ 2- - l|£ +P 7 ^ 2" - l|£ . (99) 



To calculate the first term in (I99b . we evaluate 7 when a; ^ 1 V/ 



^ ^ y-iv (100) 

(101) 



where (fTOOb follows because p^-"' + = , dToTT) follows because Ylj=i P^'"' and (fT02l) follows 
because the sum in dlOll ) is asymptotically dominated by the largest component. 

We continue to bound the first term in (l99l) 

7^2#-l|£^^p(^lp2^^i,^2^ 

Amin ^ P''^ (103) 

= (104) 

where (11031 ) is the same as (fT2l ) , hence (1104b follows. 

To calculate the second term in ( [99l ). we evaluate 7 when one or more a/ ^ 1. Consider the the two 



summations in the denominator of ( |94l ). The first one can be asymptotically evaluated as 



N 



V p - = V -J— + V pi-". 

p-(l-a„.x) |£| = N 



(105) 

max(p-i+°',pi-°") ^ p-(i-"max) 1 ^ |£| < AT 

where a' = maxaj.<iaj and a" = miUa-yiaj and dlOSI ) follows because mini p~^^"' , p^^"^") ^ p-{i-o!max)_ 
The second summation in the denominator of (|94l) can be evaluated as follows 




p-{i-"0 



(106) 



February 27, 2012 DRAFT 



24 



We now use (11051) and (1106b to bound 7 



7^ 



/OAmin |£| = 

:p3A3 i^|£|<Ar 



We thus have 



7 



Po«t 7^2« -1 



(107) 



^ P|^7 ^ 2" - 1 
< P(7 ^ 2# - 1 



S| = +P 7 ^ 2- - 1 



< 151 < TV 



P(A„,in^p-^)+F(ALn=^P"') 
P(Amin ^ P^') 



(108) 



This concludes the calculation of outage upper bound. We now proceed with the outage lower bound. 

Define the event V = {\aki\ e \/k,l} where aki is the {k,l) entry of the unitary matrix U (c.f. 
equation (|27]|). Define 

1/N 



7 



^-^7V p"-"J v^AT / pi^°i+A^ \2, 

Notice that 7 > 7 because |afc;| ^ e V A;, /. 
The outage probability is bounded as 

i 

N 



(109) 



= p[f^log(l + 7fe) ^ii) 
^ k=i ^ 

^ P r J]] log(l + 7fc) ^ i? P(T') 

^ Pr^log(l + 7) ^ i?") P(P) 
^ k=i ' 



7 ^ 2" - 1 



(110) 
(111) 
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The probability P(P) = 0(1), i.e. non-zero constant with respect to p. The proof is similar to the one 
in |[T3l Appendix A] and omitted here for brevity. We thus have 



7 ^ 2" - 1 
^ P( 7 ^ 2# - 1 
7 ^ 2" - 1 



C P(£)P 7 ^ 2« - 1 



(112) 



where (|112l) holds since P(£) = 0(1) as given by (l54l ). 

We further bound the outage probability by bounding 7 as follows. Once again consider the two 
summations in the denominator of ( 11091 ). For the first summation of ( 11091 ). we have 



N 



l — Ctj 



1 



(nl-a. + iV )2 nl~a, ^ 1^ 



,-(l-a.„ax) 



l-o„ 



|£| = iV 
1 ^ |£| < iV 



where the bound in the second line ( 11131 ) is true because 

Oj<l Qj>l Qj>l 

Using (11051 ) and (11131 ) to bound 7 Substituting back in ( 11091 ) gives: 



7 



= pAmin = iV 

pAmin 1 ^ |£| < iV 



7 



Thus the outage bound in (11121 ) can be then evaluated as we did for the upper bound 

Po«t ^P(^7^2#-1 

^ P( ^ ^ 2# - 1 



< P( 7 ^ 2" - 1 
^ ^ 2« - 1 



B\ = O^F{\B\ = 0) + 
C,0< \B\ < ivjp(|£| < N) 



(113) 



(114) 
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= P(Amin ^ p-^)O(l) +P(A^in ^ p-^)0{l) (115) 
= p-iM-N+l)_ (116) 

where dllSI l follows as a direct result of Lemma [T] From (11081 ) and (11 161 ). we conclude that the diversity 
of MIMO system using TxWF precoder and ZF equalizer is 



2) Regularized Zero Forcing Precoding: The ZF equalizer is given by (1891 ) where the composite 
channel H = HT. The received signal to noise ratio of the A;-th output symbol of the ZF filter as 



^'^ NRn{k,k) 

(117) 



l^j=l (A,+iV)2 1^1=1 \ A, ) ^"'kll 

The process of obtaining lower and upper bound has many similarities with the developments of 
Section IV-All therefore we omit many of the steps in the interest of brevity by referring to the previous 
developments. 



We begin with the outage upper bound, which is developed in a manner similar to (1961) . 

TV 



/lut = p(f^log(l + 7fc) 

^ fe=i ^ 



where 



7^2--l (118) 



p/N 



2^j=i (x^+N)^ 2^1=1 y Xi ) 



p/N 



p/N 



(119) 



< p/N 
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^ P/N 

Thus the outage in (11181) can be bounded as 

nut^pf7^2^ -1 



^log(l+7fc) 
• k=l ' 

^Prj]log(l + 7) ^i?) 

^ fe=l ' 



where 



p/iV 



Ef=i p-"^ Elli ^P^"' 

< ^4? for arbitrary 7 



p/NX 



2 

min 



(120) 



= P(Amin ^ p-°-') 

^^_i(M-7V+l)^ (121) 

We now turn to the lower bound, which is obtained in the same manner as (11121) : 

N 

-fout 

N 



7 ^ 2- - 1 (122) 



p/N 



(123) 

Let Ci = (2- - 1) eiV, C2 = CiC where ^ is a fixed positive constant (independent of p), we have 
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^ P 7 ^ 2- - 1 



A, 



Aj 



^p(^pAL„^C2)p(A,^0 



/'ALn^C^2j. (124) 

The exponential inequality (|124l) holds because P(Aj ^ ^) = 0(1), as proved in Appendix ICl We thus 
conclude: 

Remark 4: We note that the diversity of regularized zero-forcing precoder together with a zero-forcing 
equalizer can be fractional. To our knowledge this is the first instance of fractional diversity uncovered 
in the literature. 

3 ) Matched Filter Precoding: In this case, the composite channel is 

EI = HT = /3HH^. 

The noise correlation matrix is given by 

Rn = [H^H]-i = ^ [(HH^)2]-i = -1 (UA^U^)-!. 

Thus 

1 ^ 1 

i?^(A;,fe) = — V-^l^fczP (125) 
The precoder normalization factor /? = 1/ y^, where rj is given by 

N 

r/ = tr[HH^] = J]A, 

1=1 

The signal to noise ratio of the fc-th symbol of the ZF filter is 

P 



Ik 



NRn{k,k) 

p/N 



EN -v ir^Af 1 I 12 ' 
j=l ^3 2^1=1 Af I^Hl 



(126) 



Notice that the SINR 7^ in ( 11261 ) is similar to the SINR 7^ of the RZF precoding system with ZF 
equalizer given by dl 171 ). The only difference is the term A^ + A^ which, when applying the transformation 
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of Afc 



' , has no effect on the diversity analysis as detailed in the previous section. We then conclude 



that the diversity of the MIMO system applying MF precoder and ZF equalizer is the same as the diversity 
of the RZF precoder with ZF equalizer. Thus: 



d 



MFP-ZF 



l(M-iV + l). 



(127) 



B. MMSE equalizer 

The MMSE equalizer has better performance compared to ZF and is therefore widely popular. We 
investigate the diversity of MIMO systems that deploy different precoders at the transmitter and MMSE 
equalizer at the receiver. 

1) MFTx Precoding: The MFTx precoder, Tmfp , is given by ( [63] ). The MMSE equalizer for the 
precoded channel is given by 



-1 



^MMSE = 

where H = HTmfp = ^MFpHH^ and Pmfp is given by (l64l) . 
The SINR at the output of the MMSE filter is given by |[T9l 



(128) 



Ik = -^hfc 



1 



hi 



I + ^M^M 



(129) 



J kk 



where Hfc is a submatrix of H obtained by removing the k-th column, h^. 

The diversity analysis of the precoded system uses some results from the un-precoded MMSE MIMO 
equalizers |[TOl . which we quote in the following lemma. 

Lemma 3: consider a quasi-static Rayleigh fading MIMO channel H G ([^MxN ^^yj ^ ^y), the outage 
probability of the MMSE receiver satisfies 



out 



tr(I + -^H^H)-^ ^ N2- 



( ^ 1 



(130) 
(131) 
(132) 



where {A'^} are the eigenvalues of H and d'^^^^'^^ is given by ^f4\) . 
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Substituting A'^ = p we have 

1 . «'fc<l 



1 + 



(133) 



^ 11 a; > 1 



thus the term i^^^, is either zero or one at high SNR, and therefore to characterize the sum in (11311 ) 
at high SNR we count the number of ones, or equivalently the number of a'^ > 1. Hence the outage 
probability reduces to ifTOl 

Pout = ^(Y,l=\N2-^]\. (134) 

Now we apply the matched filter precoder. Similarly to (I130I ). the outage portability is given by 

Pout = P(^tr(l + ^HH^)-i ^ iV2-#^ (135) 

ErT^^^2-# (136) 

k=l^^ Nr,-^k / 

where we have used HH^ = i(HH^)^ = iUA^U^ to obtain (Il36l ).and {A^} are the eigenvalues of 
the Wishart matrix HH^. The scaling factor r/ = tr(HH^) = 

We begin with a hypothetical precoder whose transmit power is not normalized, i.e., rj = 1. The outage 
probability of this un-normalized precoder is similar to that of the MMSE receiver with no precoding at 
the transmitter, as given in (1132b . except that the eigenvalues are now squared. Thus similarly to (11331) . 
we have the exponential inequality 

1 fp2"'=-i afc<0.5 

TT^ = • (137) 

N^k \i a;->0.5 

The analysis of fTO^ then follows and we have 

d=-\\N2-Nf + {M-N)\M2~Ny\. (138) 



2 ^ 

We conclude that the un-normalized matched filter precoding with MMSE receiver results in 50% 
diversity loss compared to MMSE receiver with no transmit precoding. 

For the normalized precoder, we begin with the outage probability in (1136b . Assume ai'^ a2 - ■ ■ a^, 
the sum term in (1136b is given by 

N , N 

Y \ = V ^ 
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= y ^ , , ■ (139) 

where we have used the fact that the p^"'' is dominated by the maximum element at high SNR. 
It is easy to see that the terms of ( 11391 ) are either one or zero at high SNR, depending on whether 
p-<^N asymptotically dominates p^^^Qi qj. yj^g versa. These two cases are delineated with the threshold 
Ofc ^ 0.5max(l , + 1), or, considering that ajy is positive, ^ 0.5(aAr + 1). Thus at high SNR, 
the outage probability is evaluated by counting the ones 

N 

fc=l ^ + Nv^k 



T.— 

■ - h 



•afc>0.5(aN + l) 

J] 1 = ^) (140) 

where L = [A^2"«]. The conversion from inequality to equality in equation (11401 ) follows from 
arguments developed in ifTOl Section III-A] . 

Therefore, the outage probability is asymptotically evaluated by: 



Pout = \ (141) 
JS+ 

where P(a) is the joint distribution of the ordered ai ^ • • • ^ oat and the region of integration is defined 
as = 5 n M^+, where S is given as follows: 

• If L = A^, then we seek the probabiUty that > ^(oat + 1) for A; = 1, . . . , A^, which implies 
ttTV £ (Ij oo). Thus the integration region can be tightly represented as: 

S = {un > 1 , mill ai. > 0.5(aM + 1)| 

l<A;<Af 

• If L < A^, then we seek the joint probability that > ^{aN + 1) for = 1, . . . , L and < 
^(aAT + l) for k = L + 1, . . . ,N, implying on £ (0, 1). Thus the region of integration is represented 



as: 



5 = \aM < 1 , min ai, > 0.5(aM + 1) , max ai, < 0.5(aM + 1)) 
l<k<L L<k<N 
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Using methods similar to |[T2l and ifTOl Eq (18) - (20)], exponential equality relations can be used to 
reduce the integrand to the following: 



Pout= I Up-^''^''^''-''^^^ d{c.) 



(142) 



First we consider L = N. The probability expression is evaluated by simply taking the integral over all 
variables except on, and then taking an integral over on- 



p - / -{2N-l+M-N)aN 

J Q]v=l 

N-1 

X 

k=l 

N 



J-J ^_(2fc-l+Af-iV)(0.5+0.5a„)^(^) (^43) 
fc=l 

fc=l 

= pELi-(2fc-i+Af-7V) 

= p-*^^. (145) 



When L < N, we repeat the same integration strategy. 



p ^ ^-{2N-l+M-N)aM 
J aiv=0 

TV . 

X I ]^ _ ^-(2«-l+M-Af)(0.5+0.5a«) 

i=L+l ^ 
L 



X J]p-(2fc-l+A^-^)(0.5+0.5a„)^(^) (146) 

(27V-1+M-Af)ajv 



fc=l 

1 

P 



'ajv=0 

X ]^p-(2A.-l+A'/-^)(0.5+0.5a„)^(^) (147) 
k=l 

L 
fc=l 

.ELi-i(2fe-l+M-A') 



-\{L^+(M-N)L) 



(148) 



In deriving (1146b and (11471 ) we have used p'^^^^diat) = p~°''"' UHl. Equations (11451 ) and (1148b show 
that the system exhibits two distinct diversity behaviors based on whether L = \N2^^~\ < N. We can 
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solve to find the boundary of the two regions R = N\og j^^- To summarize: 

^MFP-MMSE _ 

'i(r7V2-#]2 + (M-iV)rM2-#l) i?>iVlog^ 

(149) 

MN otherwise 
Remark 5: The outcome is interesting for its practical implications: An MMSE receiver working with 
matched-filter precoding will suffer a significant diversity loss compared to an MMSE receiver without 
precoding, except for very low rates corresponding to R < N log 7^^, where the combination of MMSE 
receiver with matched filter precoding has exactly the same diversity as the MMSE receiver alone. 

Remark 6: Recall that R = N log is exactly the same threshold below which matched filter 
precoding (without receiver-side equalization) achieves full diversity. 

2 ) WFTx Precoding: Using the Wiener filter precoding at the receiver results in the composite channel 

e = HT = /3HH^(HH^ + p-^Niy^. 
Using the eigen decomposition HH^ = UAU^, it can be shown that 

m^m = p^v{A + p-^Niy^A^v^ (i50) 

Similar to the case of MP precoder with MMSE receiver, the outage probability of WE precoder with 
MMSE receiver is given by (c.f. (11351 )) 

Pout = Pf tr(I + ^HEI^)-i ^ N2- 

j;--^^iV2-# (151) 
k=i ^ + Thi^k / 

where {A^} are the eigenvalues of H^H and r] is the scale factor. Using ( |150b . {A^} are given by 



The scale factor rj is calculated as in (1331) 

N 



I 

Thus the outage probability can be written as 

N 



Pout = r(f2^k>N2-A (153) 
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where 



1 P P 



Ik 



A 



where we define Vk = We now proceed to express both p~^r] and Vk in terms of {a^}, the 

exponential orders of {A^}. 



1=1 



{p-^Xi + Ny ^(pi-". + iv)2 



cei>l ai<l 

observe that all the terms in (1154b have negative exponent. Using (I152I ). 



Vk 



1 

1 p2(l-a.) 



1 afe < 1 

Ofc > 1 

From (1154b and (11551 ). we see that when < 1 then i;^ + p^^r] = = 1. On the other hand, when 



(155) 



afc > 1 then 



«1>1 «!<! 



o:i>l ai<l 
l^k If^k 



= p-\ (157) 
where (11561) follows because > 1. Thus we have 

Ik = 1 / = <^ (158) 

and p^^T? has negative exponent thus vanishes at high SNR. 

Observe that (11581 ) is similar to (1133b which corresponds to the case of the MMSE-only system (i.e. 
with no precoding). Thus substituting ( 11581 ) in the outage probability ( |153b and repeating the same analysis 
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of the MMSE-only system as in llTOl . we conclude that the diversity of the MMSE receiver when using 
WFTx precoding is the same as the diversity of the MMSE receiver with no linear precoding, which is 
given by ( [74l ). 

3) RZF Precoding: Using the Regularized Zero Forcing precoding at the receiver results in the 
composite channel 

e = HT = /3HH^(HH^ + cl)'^. 



where c is a fixed constant, (3 = l/rj and t] is given by (|33 

N , N 



k ^ ^ 

Similar to (11511 ). the outage probability of RZF precoder with MMSE receiver is given by 

N 



and 



Ik 



where {A^} are the eigenvalues of H^H given by 



\2 -2ak 

Afc = 7^^^ = T^^— k=l,...,N (160) 



(Xk + C)2 (p-"^ + c)2 

Notice that at high SNR we have 



^ -a, 



1=1 

^k 



C2 



Thus the SINR is given by (c.f. (fT39l )) 



ik 



k = l,...,N 

which are the same terms as in (11391 ). implying that the outage probability of the MMSE receiver working 
with the regularized zero-forcing precoder is asymptotically the same as the outage probability of the 
MMSE receiver working with the matched filter precoder. This means: 

^RZFP-MMSE _ ^MFP-MMSE 
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5 10 15 20 25 30 



SNR 



Fig. 5. MF and regularized ZF preceded 2x2 MIMO system for rates (left to right): R — 1.9, 2.5, and 4 b/s/Hz. 




SNR 



Fig. 6. MIMO system with matched filtering preceding and ZF equalization for rates (left to right): i? = 1, 2, and 4 b/s/Hz. 



VI. Simulation Results 

This section produces numerical results for the outage probabilities of ZF, regularized ZF (RZF), 
matched filter (MF) and Wiener precoding systems. Figure [3] shows the outage probabilities of the ZF 
and Wiener-filter precoded 2x2 MIMO systems. The diversity in the case of the ZF case is the same 
as the one predicted by the DMT. In the case of Wiener precoding, the diversity is the same as the 
one predicted by the DMT for high rate (R) values and it departs from the DMT for low rate values. 
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SNR 



Fig. 8. 2X2 MIMO system with Wiener filtering preceding and MMSE equalization for rates (left to right): R = 1.5, 3, and 
4 b/s/Hz. 



A complete diversity characterization is given by (1741 ) which is similar to that of the MMSE MIMO 
equalizer |10]. Figure |4] shows outage probabilities for a 3 x 3 MIMO system with Wiener precoding. 
The diversity for the rates R = 1.5, 4, and 5 b/s/Hz is 9, 4 and 1 respectively. Figure |5] shows an error floor 
for the regularized ZF and matched filtering precoded 2x2 system at high rates. However we observe 
that the maximum diversity is achieved for any rate R < 2 (c.f. Equation (l58l)). Figure [6] shows outage 
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SNR 



Fig. 9. 3X3 MIMO system with Wiener filtering preceding and MMSE equalization for rates (left to right): R = 1.5, 4, and 
5 b/s/Hz. 



0) 
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— 3 b/s/Hz 




' 2:5 bSl1r^'^~^-^---_J_^ 
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;':>:)SH/ 


1 5 Ws/hX. \ 




MF precoder/MMSE equalizer 




— e — MMSE equalizer only 




1 ^ _ ^ ^ 





5 10 15 20 25 30 

SNR 

Fig. 10. 2X2 MIMO system with MF preceding and MMSE equalization system for rates (left to right): R — 1.5, 2.5, and 3 
b/s/Hz. 

probabilities for a 2 x 2 and a 3 x 3 MIMO system with matched filter precoding and ZF equalization. 
The observed diversity values are consistent with Eq. (11271 ). Figure |7] shows outage probabilities for a 
2x2 and a 3 x 3 MIMO system with Wiener filter precoding and ZF equalization. Figure [8] and Figure |9] 
show outage probabilities for a 2 x 2 and a 3 x 3 MIMO system, respectively, with Wiener filter precoding 
and MMSE equalization. The diversity for the 3x3 system is the same as the diversity of the Wiener 
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filtering precoding-only (c.f. Figure |4]). 

Figure [10] shows the outage probability of a 2 x 2 MIMO system with matched filter precoding and 
MMSE equalization, which is consistent with Eq. ( |149b . We also plot the outage probability of the MMSE 
MIMO equalizer (without any precoding) for comparison. 

VII. Conclusion 

Linear precoders provide a simple and efficient processing, and have been shown to be optimal in 
some scenarios fSj, Q, Q. This paper studies the high-SNR performance of linear precoders. It is 
shown that the zero-forcing precoder under two common design approaches, maximizing the throughput 
and minimizing the transmit power, achieves the same DMT as that of MIMO systems with ZF equalizer. 
When a regularized ZF (RZF) precoder (for a fixed regularization term that is independent of the signal- 
to-noise ratio) or matched filter (MF) precoder is used, we have d{r) = for all r, implying an error 
floor under all conditions. It is also shown that in the fixed rate regime RZF and MF precoding achieve 
full diversity up to a certain spectral efficiency, while at higher spectral efficiencies they produce an error 
floor. If the regularization parameter in the RZF is optimized in the MMSE sense, the RZF precoded 
MIMO system exhibits a complex rate-dependent behavior. In particular, the diversity of this system (also 
known as Wiener filter precoding) is characterized by d{R) = [A^2~ n]^ -|- (m — N) \N2^^] where M 
and are the number of transmit and receive antennas. This is the same behavior observed in linear 
MMSE MIMO receivers ifTOl . Various results for the diversity in the presence of both precoding and 
equalization have also been obtained. 

Appendix 

A. Pairwise error Probability (PEP) Analysis 

In this section we perform PEP analysis for the the zero-forcing (ZF) and the regularized ZF (RZF) 
precoding systems. The presented analysis can be easily extended to all other precoding systems. The 
basic strategy is to show the SNR exponent of outage probability bounds the SNR exponent of PEP from 
both sides The PEP analysis follows from |[T4i . |[TOl . with careful attention to the system model given 
by Equation ([1). 

The lower bound immediately follows from lfT4l Lemma 3] by recognizing that although it was 
developed for SISO block equalization, nowhere in its development does it depend on the number of 
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receive antennas, therefore we can directly use it for our purposes: 

Perr ^ Pout- (162) 

The upper bound on PEP for the ZFIKLV precoding systems receiver is developed using the union 
bound. Denote the channel outage event by O and the eiTor event by E. The PEP is given by 

P,n = P{E\0)Pout + P{E,0) 

^Pout + P{E,0). (163) 

In order to show that Pout dominates the right hand side of (11631 ). it is shown in ifTOll that the probability 
P{E, O) can be bounded as follows using the union bound 

p/N 

F{E,0) ^ 2™e ^ ^ p-^^^ (164) 

where / is the codeword length and (T?(fc) is the variance of the interference plus noise signal h in the 
k-th receive stream The proof of |[T4ll does not depend on the codeword length for both upper and lower 
PEP bounds. The bound are tight and were confirmed by simulations for outage and error probabilities. 

We now show that a similar proof holds for regularized zero-forcing (RZFP). Recall that the outage 
probability of the RZFP can be upper bounded by (|6TI ) 

Pout ^ ^ e) ^ (165) 

We will use Pj^^ to further bound ( 11631 ). Moreover P{E, O) can be upper bounded by bounding the 
noise variance cr'^{k) in (1164b 

al{k) = Pi + Pn<PT + l (166) 

where we have used the noise power P„ = 1, and bound the interference power by the total received 
power Px- We will first consider the case of RZF precoding since the case of ZF precoding can be easily 
deduced from RZF by substituting setting the regularization parameter c = 0. For the RZF precoding 
system we use the Pt given by (l34l) which can be simplified in a way similar to earlier sections 

P 

* 1141 analyzes linear receivers so n is the fc-th output filtered interference plus noise signals. By symmetry assumption all 
the equalizer outputs have equal noise variance. 
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E 



1 p ^ p 



Ipi-"™. (167) 



Using the union bound ( 11641 ). 



P(i?,0) ^ <; ^ (168) 



Since the exponential function dominates polynomials we have 

lim 

and 



which in turns gives 

P{E, O) ^ p-^'^ . (169) 
Using ( 11651 ) and ( 11691 ). the PEP given by ( 11631 ) is bounded as 

P,n ^ Pout + P{E, O) 

= 

out 

= p-'^°"*. (170) 

therefore d ^ dout which concludes the proof for the RZF system. 

For the ZF precoding system, it can be directly shown that a similar proof holds for both ZF precoding 
designs. 
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B. Proof of Eq. ^ 
Recall that 



N 



* 12 ■ 

k=2 



All terms of ih the common factor , — Thus we have 

1 



'^b= (r^ + T^ + r^ + T^ + "' + r^)- ^^^^^ 

Observe that all the terms of tpb are distinct except for the first two. 
We now bound the probability Fi^ip ^ P~)- 

P(V' ^ p^) ^ IP(V' ^ I < c) P(^A < c) 

^ P(c ^ < C) 

= P(V' < c) (172) 
Using Tp = Tpatpb we can further bound (11721 ) 

P(V; < c) = P(VaV'fe < c) 

^ P(V'„V6 ^ C |Va < Cs) P(^a < C2) 
^P(C2V'6^C) P(^a<C2). 

We thus have 

P{i'i^p-)^^{A^ c') ni^a < C2) (173) 

and c' = c/c2- 

We now evaluate the two probabilities in the right hand side of ( I1731 ). The first probability P(V'6 ^ 
c') = 0(1). The proof easily follows from [13, Appendix A] with the observation that this proof holds 
even when the two first elements of ipi, are the same. The second probability F{ipa < C2) is evaluated as 
follows. Let q = lui^'p. We use the following distributions from [9, Appendix A] 

f(q) = (N-l)il-qf-^ O^q^l 
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then 



b < C2) = FiQ > -) 
fiq) dq 



C2 



N-2 



(174) 



Observing that ( 11741) is not a function of p concludes the proof. 



C. Proof of F{Xi ^ = 0(1) for any I 

Define a Wishart matrix W using the Gaussian matrix H. 

H M > N 
H^H N i^N 



W 



Let n = max(M, A'') and m = min{M,N). The matrix W is m x m random non-negative definite 
that has real, non-negative eigenvalues with Ai ^ • • • ^ AmO. The joint density of the ordered eigenvalues 



is im 



/(A) = K^^e- ^- n ^"""^ n(^^ - 

i i<j 



(175) 



Thus the marginal distribution of A; is given by ifTTI 

/a,(AO = J ...J /(A) dX2...dX. 

m 

-^^^.(AO^Ar-e-^^ 

where 

V'fc+ilA) 



i=l 



{k + n — m)\ 



1/2 



L^-'"(A), A; = 0,...,m-1 



where L^-*" 



„\ 2_ X m—n _d^{ —X n—m+k\ 

order k. 

We now compute f{Xi ^ 



(with Lo = 1) is the associated Laguerre polynomial of 



i=l 



1 



1^ ni 
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1 



m{n — m)l 
1 



m{n — my. 



n—m N 

e"^' ^n{n-l)...{n-k+ J 
fc=i ^ 



(cn-m 



m{n — m)\ 

n—m 



fc=i 



(176) 



(177) 



where (11761 ) follows from ||22l Section 2.32]. The right hand side of Equation (11771 ) is a non-zero constant 
bounded away from zero. This concludes the proof. 
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